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I. INTRODUCTION 

In this report, we extend to cloud physics the work done in (5-9) for single and 
multiple scattering of electromagnetic waves. We consider the scattering of light (visible 
or i nfr ared) by a spherical cloud represented by a statistically homogeneous ensemble 
of configurations of N identical spherical water droplets whose centers are uniformly 
distributed in its volume V. The ensemble is specified as in (8), by the average number 
p of scatterers in unit volume, and by p/(R) with /( R) as the distribution function 
for separation R of pairs. The incident light, <f> 0 = aoe lfc °* r , a plane electromagnetic 
wave with harmonic time dependence, is from outside the cloud. The propagation 
parameter k Q and the index of refraction rj 0 determine physically the medium outside 
the distribution of scatterers. 

We solve the interior problem separately to obtain the bulk parameters for the 
scatterer equivalent to the ensemble of spherical droplets (2—5). With the interior 
solution or the equivalent medium approach, the multiple scattering problem is reduced 
to that of an equivalent single scatterer excited from outside illumination. A dispersion 
relation which determines the bulk propagation parameter K and the bulk index of 
refraction 77 of the cloud is given in terms of the vector equivalent scattering amplitude 
Q and the dyadic scattering amplitude g of the single object in isolation. 

Based on this transfer model we will have the ability to consider clouds composed 
of inhomogeneous distribution of water and/or ice particles and we will able to take 
into account particle size distributions within the cloud. We will also be able to study 
the effects of cloud composition (i.e., particle shape, size, composition, orientation, 
location) on the polarization of the single or the multiple scattered waves. Finally, this 
study will provide a new starting point for studying the problem of lightning radiative 
transfer (3-4). 

In general, we work in spherical coordinates. We use bold face or an arrow to denote 
a vector or a vector operator. A circumflex indicates a vector of unit magnitude. A 
tilde on the top of a letter denotes a dyadic (second rank tensor). For brevity, we use 
[5:4] for equation 4 of Ref. ( 5) etc. 

II. MATHEMATICAL MODELING/SOLUTION INSIDE THE CLOUD 
FOR OUTSIDE INCIDENCE 

The solution inside the cloud for outside illumination corresponds to the multiple 
scattering of a plane electromagnetic wave by an ensemble of configurations of N iden- 
tical spherical water droplets. To obtain the solution inside the cloud, we consider first 
the single scatterer in isolation, second a fixed configuration of N identical scatterers, 
and third an ensemble of the above-mentioned configurations. 

For an incident plane electromagnetic wave <p = aie ,#Cl ‘ r , = k 0 rj ly and rj 1 being the 
complex relative index of refraction for the host medium inside the cloud but outside 
each droplet, the total outside solution for the single scatterer in isolation (outside 
the single water droplet but still inside the host medium) tp = <p + u satisfies the 
following differential equation obtained from Maxwell’s equations after suppressing the 
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harmonic time dependence factor e ,u,i 

|v x V x J ip = 0, V ■ ip = 0. 


[ 1 ] 


The solution inside the single spherical water droplet in isolation xpi n satisfies [lj 
with replaced by « a . Here, /c 3 = r] 3 = h with t? 2 being the complex relative 

index of refraction for the medium inside the spherical water droplet. The propagation 
parameters and k 3 correspond (within the distribution of identical spheres) to the 
media outside and inside the water droplet respectively. 

Similar to Twersky (7), we have 

ip = a a e , * l ' r + - r'|),u(r')},u(r) = {/», uj = 

/ [(* X a) • (v x u ) - (v x h) • (n x u)] dS(r'). P1 

Here, h — — r^l), h(x) = and I being the identity dyadic. It 

is important to note that r, and r' denote the observation point and a point on the 
surface S or in the volume v of the water droplet respectively. 

Asymptotically (acjT >> i) we can write 

u( r ) = *(« 1 r)g(r,« 1 : aj, g(?)= I. • g(r). [3] 

Here, I* = (i - rr) is the transverse identity dyadic and a ! • k x = 0. The spectral 
representation of u is 

u ( r ) = ~ J e '* lc t g(f )df2(6 c , tp e ) y r > (r • r'), k 1c = K ac r c (0 c , <p c ), [4] 

C 

and the single scattering amplitude u(r') j can also be 

evaluated from Mie scattering theory. 

Now, we consider a fixed configuration of N identical scatterers with centers located 
by r m(m=i,2,3,...,A r )- The total outside field 

N 

*( r ) = <fc*) + ^2 U m (r - r m ), U m (r- r m ) ~ ft(K,|r- r m |)G m , |r-r OT |-»oo. [5] 


Equivalently, for the scatterer located at r*, we use the self-consistent approach of 
(6—7) to obtain the total outside configurational field 

*<(>•) = ^(r) + 53' u ro ( r - r m ) + u ,( r - r,), = £ . [6] 
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Using [6] and the general reciprocity relation 4, } = 0 for any arbitrary direction 
of incidence and polarization a x , we derive as in (2) the self-consistent integral equation 
for the multiple configurational scattering amplitude 

G 1 (r) = g 1 (f,*,)-a.e iiir, + £'/ g,(f,f e )'G m (f e )e““' R “, [T] 

c 

with Rt m = r t - r m , / = £ f dn c , g(r, k ,) • a, = g(f, k, : aj, and the magnitude of 

the separation distance |Rtml I s bounded above by the diameter D of the cloud. 

We take the ensemble average of [7], use the quasi-crystalline approximation of Lax 
(2), the equivalent medium approach and Green’s theorems, to obtain (7) the dispersion 
relation determining the bulk parameters 

k) = 0]}+, / [/(R)-i]e-' RR 0«R, [8] 

Vv-V 

where dR denotes volume integration over (Vfo — v). Here, Q is the equivalent scattering 
amplitude and U is a radiative function defined by U = f g(r, r c )*C? f^icl Kj e** ic and 

Co = Kx/Airi. The bulk propagation parameter K = K x r} with t) being the bulk index of 
refraction, and {[/,$]} = / [fd n g - gd n f]dS is the Green surface operator with outward 

5 

unit normal from v. Equation (16) solves formally the interior problem for the cloud 
with outside illumination. 

III. BULK PARAMETERS AND LEADING TERM APPROXIMATIONS 

To simplify [8], we force the model to neglect all phase transition effects (1), and 
to take only into account pair interaction due to central forces. If the inter-droplet 
potential is negligeable, we can choose /( R) to be always equal to unity. Hence, [8] 
is reduced to 

(K a - *5)1, + (£)g(f,K)]-£(*,| k) = 0. [9] 

In [9], let r = K. In addition, because optical scattering from a cloud is highly 
forward peaked (7), we neglect back scattering and reduce [9] to 

(K a -*5 )I.(k) + (£)g(K,K)] e(K| k) =0. [10] 

If the scatterers preserve the incident polarization (7:68), we have from [10] 
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Equation [11] determines the bulk propagation parameter K and the bulk index of 
refraction rj of the equivalent medium for the bounded distribution of the spherical 
water droplets. 

IV. CONCLUSION 

The multiple scattering problem has been reduced to that of a single equivalent 
scatterer in isolation. Formulae are given for the bulk propagation par am eter K and 
the bulk index of refraction 77 of the equivalent medium. The results axe quite general 
in nature and can be extended to non-spherical geometries. Also, they can be applied 
immediately to the problem of pulsating optical point sources arbitrarily distributed 
throughout a scattering medium. When /( R) — 1/0, [8] can be approximated or 
solved numerically. 
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